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• Abstract 

^ I Since the Standard Model most probably cannot explain the large value of CP 

asymmetries recently observed in D-meson decays we propose the fourth quark-lepton 
generation explanation of it. As a byproduct weakly mixed leptons of the fourth 
generation make it possible to save the baryon number of the Universe from erasure 
by sphalerons. An impact of the 4th generation on BBN is briefly discussed. 



1 Introduction 

Recently LHCb collaboration has measured the unexpectedly large CP violating asymmetries 
in D — 7- TT+TT^ and D — )■ K^K~ decays [1]: 

AA^f = Acp{K+K-) - Acpin^TT-) = [-0.82 ± 0.21(stat.) ± 0.11(syst.)]% , (1) 

where 

^^^^^ ^ ) = nD^-^n^n-) + TiDO^n^n-) 

and Acp{K^K~) is defined analogously. 

This result was later confirmed by CDF collaboration, which obtained [2]: 

AAg^^ = [-0.62 ± 0.21(stat.) ± 0.10(syst.)]% . (3) 
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Figure 1: Quark diagrams describing D — > tt+tt^ decay in the Standard Model. A wavy 
line denotes VT-boson, a curly line - gluon. 



The most important question concerning experimental results ([T]) and (13]) is whether in 
the Standard Model the CP-violation (CPV) in these decays can be as large as 0.5% - 1%. 

In the Standard Model the CPV in D{D) — )■ vr+vr" decays originates from the interference 
of the tree and penguin diagrams shown in Fig. 1. For D{D) — )■ K^K~ decays rf-quarks in 
these diagrams should be substituted by s-quarks. 

It is convenient to present the penguin diagram contribution to D ^ tt+tt^ decay ampli- 
tude in the following form [3]: 

= VcdV:Af{m,) - /K)] + v,,v:,[f{m,) - /K)] , (4) 

attributing the first term to the tree amplitude and considering the second term only as the 
penguin amplitude. 

In the case of — )■ K^K^ decay the following presentation is useful [3]: 

= KsCl/K) - f{^d)] + V,,V:,[f{m,) - /(m,)] , (5) 

where the first term is attributed to the tree amplitude while the second one is the penguin 
amplitude. 

Denoting the absolute values oi D ^ tt+tt^ decay amplitudes by T and P we get: 



1 + ^Jis-i) 
T 

T 



(6) 



where S stands for the difference of the strong interaction phases of the tree and the penguin 
amplitudes, while 7 ~ 70° is the phase of Vub (the product VcdV*^ as well as Vcb are practically 
real in the standard parametrization of the CKM matrix). 
From eq. (E]) for the CPV asymmetry we obtain: 

P 

Acpin^n ) = 2— sin5sin7 , (7) 

where in the denominator of we neglect the terms of the order of P/T and [P/T^ which 
is a very good approximation because P/T ^ \VcbV*fj\/Vcd ^ 1- Here sin 7 is close to unity 
and we use this value in what follows. 
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Let us present an argument demonstrating that 6 can also be close to 90*^. The tree 
diagram gives dominant contribution to the D — > vrvr decay rates. The corresponding to it 4- 
fermion Hamiltonian has parts with isospin 1/2 and 3/2. That is why the produced vr-meson 
may have isospin zero or two. So three decay probabilities, — > tt+tt^, — ?■ tt+tt", and 

— )■ 7r°7r°, depend on the absolute values of the decay amplitudes Aq and A2 and their 
strong phases difference 60 — 62- From the experimentally measured branching ratios [1]: 

Br(L)+ ^ tt+tt") = [12.6 ± 0.9] ■ 10"^ , Br(D° ^ 7r°7r°) = [8.0 ± 0.8] ■ 10"^ , 

Br(D° ^ TT+TT-) = [13.97 ± 0.26] ■ 10"^ (8) 
we find for the phase difference of the amplitudes with / = and 1 = 2: 

\6o-62\ =86° ±4° . (9) 

In eq. ([7]) 6 stands for the difference of the strong phases of penguin amplitude which has 
1 = 1/2 and produces pions with 1 = and tree amplitude, which has parts with J = 1/2 
and 1 = 3/2 and produces pions with 1 = and 1 = 2, that is why 6 6q — 62- Nevertheless 
eq. ([9]) demonstrates that 6 can be large, and so we substitute sin 5 = 1 into eq. (171). 

In the limit of ?7-spin (d -H- s interchange) symmetry the tree amplitude of D{D) — )• 
K^K~ decay differs by sign from that of D{D) — )■ tt+tt^ decay, while the penguin amplitudes 
of these decays are equal, that is why 

Acp{K+K-) = -Acp{7r+n-) . (10) 

However since [1] 

Br(D° ^ K+K~) = [39.4 ± 0.7] ■ 10"^ , (11) 

we obtain from eq. ([8]) that \Ax+k- /At^+t^-I 1.7 and f/-spin symmetry is heavily broken 
in D decays. Nevertheless let us suppose that ffTOj) is not badly violated, so finally we get: 

AAcP = A^. (12) 
Now let us try to understand if in the Standard Model we can obtain 

^ = 1.8■10-^ (13) 
which is needed to reproduce the average value of the LHCb and CDF results. 



2 Z) — tttt: charmed penguin 

Though the four-fermion quark Hamiltonian responsible for these decays is known, strong 
interactions does not allow to make an exact calculation of the decay amplitudes. What can 
be done is an estimate of the decay amplitudes with the help of factorization. Let us start 
from the tree diagram shown in Fig. la which dominates in the decay amplitude: 

T = ^V;d< vr+TT- 1^7,(1 + 75)cM7a(l + 75)^^1^° > X 

X |^[a,(m,)/a,(MH.)]-^/^ + i[a,(m,)/a,(MH.)]^/'} , (14) 
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where the last factor originates from the summation of the gluon exchanges in the leading 
logarithmic approximation. Substituting into it 6 = 11 — 2/3Nf = 23/3, as{Mw) = 0.12, 
tts("^c) = 0.3 we find that the factor in the curly brackets is close to one, {...} = 1.1. 
Factorizing the decay amplitude we obtain: 

T = l.l^V;d<7r+|M7„(l + 75)c/|0><7r-|(i7,(l + 75)c|D°>= 

= ll^KdAA;i„[/^(0)(p + fc2)a + /:(0)(p-A;2)„] = l.l^Kd/./T(0)m^, (15) 

where ki and k2 are the momenta of the produced 7r-mesons, p is the D-meson momentum 
and we neglect in comparison with mjj. 

The value of the D° — )■ vr+e"'"!/ transition formfactor at = can be found in ref. [1]: 

/T(0)|Kd| =0.152 ±0.005, /^(0)=0.66, (16) 

and for the decay width we obtain: 

ptheor. _ Gp (l.lK:d/+(0)/,r^D)^ _ a o 1 n9o-l ^ 7^ 

- ~2 16^n - ■ 10 s , (17) 

where f-,, = 130 MeV was used. 

From the branching ratio of the — > tt~^7t~ decay ([8]) and D°-meson mean life, t^o = 
0.41 ■ 10-12 gn^. 

r^^^...- = 3.4 ■ lO^s-i , (18) 



So the naive factorization overestimates the decay amplitude by the factor y 6.2/3.4 ^ 1.4. 
Calculating the D — )■ K^K~ decay probability we obtain: 



ptheor 
D^K+K- 



r*^h-^^_ = 12.2 . lO^s-i , (19) 



where we substituted fxl f-n = 1-27 and /:f^(0) = 0.73 taken from ref. [1]. 
From eq. (fTTl) it follows: 

^""Zk^k- = 9.6 • lO^s-i , (20) 



so the factorization overestimates the decay amplitude by the factor y 12.2/9.6 = 1.1. 

We see that in the case of the tree diagrams the accuracy of the factorization approxi- 
mation is very good. 

Let us make a brief remark on the D — )■ K^K^ decay. At the tree level it proceeds 
through the diagram with VT-boson exchange in t-channel, so it should be suppressed. Even 
more, cu — )• dd and cu — )■ ss amplitudes interfere destructively and in the [/-spin symmetry 
limit their sum is zero [5]. According to experimental data |?]: 

Br(L)° ^ K^K^) = 4Br(L)° ^ 2K^s) = (6-8 ± 1.2) . 10"^ , (21) 

which is approximately 6 times smaller than Bt{D — )■ K^K^). It means that the decay 
amplitude is smaller than that to charged kaons by factor 2.5. This unexpectedly small sup- 
pression may indicate that large distance effects like — K*^K*~ — )■ K^K'^ rescattering 
can be important. 
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The four-fermion QCD penguin amplitude which describes D — )■ tt+tt decay looks like: 



H{P) 



Gp 

72 1271 



asirric) , rrib 

in 

127r 



(M7a(l + 75)Ac)(rf7aArf) = 
In [(m7,(1 + ^^)d){d^^{l + 75)c) 



- 2m(1- 75)^/^(1 + 75)c]- , 
where A are the Gell-Mann SU(3) matrices and we use the Fierz identities: 



(22) 



AafeAcd 

^7a(l + 75)</'X7a(l + 75)^ 

V'7a(l + 75)<^X7a(l-75)^ 



^7a(l + 75)^X7^(1 + 1^)^ : 

-2?/'(l - 75)r/x(l + 75)<^ • 



Also the identity < ■K^\uaOdb\Q >= l/36ab < n^\uOd\0 >, where O = 'jalb or 75, was used. 

Calculating the matrix element in the factorization approximation with the help of the 
equations of motion for quark fields we find: 



V2 



127r 



In 



rrLrJ 9 



1 + 



2m2 



(23) 



Dividing it by the experimental value of the tree amplitude and using Eq.f lTSj) we obtain: 



P/T = - — 



\Vctv:, 








12n 



rrir 



1 + 



2mt 



Substituting \Vcd\ = 0.23, \Vub\ = 3.9 ■ 10"^ K 



cb 



41 ■ 10-3. 



a.Arrir 



rric = 1.3 GeV, + nid = 6 MeV we come to: 

P/T ^ 9 ■ 10"^ 



(24) 

0.3, nib = 4.5 GeV, 
(25) 



Comparing it with eq. f lT^ we see that in order to fit the experimental data on AAqp the 
penguin amplitude should be enhanced by the factor 20 in comparison with what factoriza- 
tion gives. Concerning the tree amplitudes, we have found in this section that factorization 
result differs from the experimental value by the factor 1.4 in the case of D ^ tt^ti~ decay 
and by 1.1 in the case of D — )■ K^K~ decay. In the next two sections we will study how 
accurate is the factorization approximation to the penguin amplitudes in B- and ii'-meson 
decays. 



3 B ^ ttK: beautiful penguin 

Bu — 7i~^K^ decay is described by the penguin amplitude shown in Fig. 2. 
The Hamiltonian responsible for this decay looks like: 

H = ^VtbV*[csOs + c^O, + c,05 + cgOe] , (26) 
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Figure 2: i?„ — t- ti^K^ decay proceeds through the penguin amphtude only. 



VtbV-i*g is substituted for VcbV*^ + VubV*^ (the contribution of a loop with the virtual t-quark 
is negligible) and 

O3 = s7o(l + 75)M7o(l + 75)^ 
O5 = S7„(l + 75)6ci7„(l - 75)(i 

Oq = Sa7a(l + 75)Mc7a(l - 75)'^a , (27) 

where a,c = 1,2,3 are the color indexes. 

Using the Fierz identities as well as < K^\saOdb\0 >= ^6ab < K'^lsOdlO > identity we 
obtain: 

H = ^VtbV*[a,s^ail + 75)dd^ail + 75)fe - 2a6s(l - ^5)dd{l + 75)6] , (28) 

where 04 = |c3 + C4, = |c5 + Cg. Calculating the matrix element in the factorization 
approximation we obtain: 



M = ^VtbV;jKU{0)ml 



CI4 -|- Clg 



2m^ 



(29) 



mbms_ 

where in the leading logarithmic approximation the following approximate equation is valid: 



aA = a^ = T^ln ~ -0.03 , (30) 

127r V rrib 



while at NLO approximation from Table 1 of [6j we obtain: 04 = —0.031, ag = —0.042. 
Substituting = 100 MeV, nib = 4.5 GeV we find: 

r;t^..xo) = ^|^/i/^(0)m|(0.076)^ = 4.1 ■ lO^s"^ , (31) 
where V^s = 39 • 10^^ and /+(0) = 0.25 from [4J was used. The experimental result is: 

rS-...^o) = 14 ■ lO^s-^ , (32) 
So, the factorization result is enhanced by the factor 



P/Pfact = V14/4.1 = 1.8 . (33) 

The numerical value of the penguin amplitude is important in the calculation of CP 
asymmetries in — )■ irK and i? — )■ vrvr decays [Tj. 
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4 K ^ tttt: strange penguin 



s — 7- (i penguin transition changes the isospin by 1/2 in this way explaining the famous 
A/ = 1/2 rule m K ^ nn decays. The calculation of the Ks — )■ 7r~^7r~ decay amplitude 
generated by the penguin transition using the factorization underestimates the amplitude 
by the factor 2-3 according to refs. [6], |8] . 

In view of the results for B and K decays we can cautiously assume that for D — > n^n^ 
decay the factorization calculation underestimates the penguin amplitude at most by factor 
5 leading to: 

(AA*^^r)^^,^0.2% . (34) 

Thus the following alternative emerges: either the experimental results are wrong or New 
Physics is found. Of course we cannot determine what kind of new particles and interactions 
are responsible for large CPV asymmetry in D — )■ tt+tt^ {K~^K~) decays. However, in the 
next section we will propose the straightforward generalization of the Standard Model in 
which large CPV in D decays can be explained. 



5 The fourth generation: enhancement of CPV in D 
decays 

As it was stated in paper [0] the introduction of the fourth quark-lepton generation may 
easily remove Standard Model upper bound ( IMl) matching the experimental results [H E]. 
In the case of the fourth generation the additional term with the intermediate b' quark 
should be added to the expression for the penguin amplitude. In this way expression (jl]) is 
substituted by: 

= VcdV:,[f{md) - /K)] + KfeKJ/K) - /K)] + V,f,'Vub'[fim,,) - /(m,)] , (35) 

where the unitarity of 4x4 quark mixing matrix is used. According to the experimental 
constraints from the direct searches of the fourth generation quarks b' should weigh several 
hundreds GeV, that is why fimy) is small and can be neglected just as it is done with 
t-quark contribution to 6 — )■ s penguin, see the remark after Eq. (1261) . In order to enhance 
SM contribution to the penguin amplitude we should suppose that the term Kfe'Kife'/("^s) 
dominates. 

Then the enhancement of Acp in the case of the fourth generation is equal to: 



P4 




\Vct'V:,,\ 


sin(argKfe/1/V) 


PSM 




\VcbVub\ 


sin 7 



3 • 10-4 , , 

where in the last equality we use the allowed values of the product iKfe'V^Vl sin(arg Vc^'Kife') 
taken from Fig. 1 of paper [in] 0- So we see that the enhancement necessary to describe the 
experimental data on AAcp can be achieved in the case of the fourth generation. 



^Let us stress that the logarithmic {log (mw / "m-c)) enhancement originates not from the diagram with 
the intermediate b' quark but from the term f{ms). 
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6 Saving baryon number by long-lived fourth genera- 
tion neutrino 



If weakly mixed particles exist, then the sphaleron processes can create the baryon asymme- 
try of the universe [11]. As it is noted in ref. [12], the long-lived fourth generation particles 
save baryon asymmetry generated in the early universe from erasure by the sphaleron tran- 
sitions. The sphaleron transitions conserve B — L, thus, if in the early universe Bq = Lq ^ 
is generated, then the final baryon and lepton asymmetries being proportional to B — L are 
completely erased. If the fourth generation particles weakly mix with three quark-lepton 
generations of the Standard Model, then two additional quantities are conserved: i?4 — L4 
and L — 3L4, where B4 and L4 are the densities of baryons and leptons of the fourth gen- 
eration, while B and L are the densities of baryons and leptons of three light generations. 
In ref. pr2] initial asymmetries Bq = Lq = 3A and B^ = L\ = were chosen and since 
L — 3L4 = 3A 7^ 0, the total baryonic number density, B + B4^, being proportional to a linear 
superposition of conserved quantities is nonzero at the sphaleron freeze-out temperature. 
After the sphaleron freeze-out B + B4 is conserved in comoving volume and is equal to the 
present day baryon density of the Universe. However, if heavy baryons of the 4th generation 
do not decay prior to big bang nucleosynthesis (BBN), the light baryon number density at 
BBN could be different from that determined from the angular fluctuations of CMB. The 
impact of this effect on the light element abundances is discussed below. 

For such a scenario to occur the lifetimes of the fourth generation quarks and leptons 
should be at least larger than the universe age at the sphaleron freeze-out: r4 > Mpi/Tgpj^ ~ 
10~^° s. For the mixing angles in the case of b' — )■ (c, u)W decay it gives 6 < 10^^ [12j, much 
smaller than what we need to explain the large CPV in D-decays, see eq. (l36i) . 

So in our case quarks of the fourth generation should be much stronger mixed with quarks 
of three light generations. However, let us suppose that leptons of the fourth generation 
are weakly mixed with the leptons of three light generations. Let us introduce the total 
baryon density, B' = B + B4, and take the initial conditions analogous to those in ref. |12] : 
Bq = Lq = 3 A and L4 = 0. We can choose four independent chemical potentials as: fj,w, 
= + ^u^, + yUj/r) which are the chemical potentials for the upper type quarks, 
H^-bosons, 4G neutrino and sum over all SM neutrino chemical potentials (see Appendix). In 
the limit /ij/T <C 1 the baryon and lepton densities are linear combinations of these chemical 
potentials with the coefficients which depend on the ratio of masses of the corresponding 
particles to the temperature. We will take into account the masses of ly-boson, t-quark, 
t'- and 6'-quarks of the fourth generation and the fourth generation leptons and the 
masses of all the other components of the primeval plasma can be neglected in comparison 
with Tsph. 

Finally we have four equations for four unknown chemical potentials: two quantities are 
conserved under the sphaleron transitions; we can choose them as 

B'-L~L4 = , 

L-3L4 = 3A . (37) 

The third equation is that of the electric neutrality of the primeval plasma, Q = 0, and, 
finally, the sum of the chemical potentials of all the particles which are converted into 
nothing by sphaleron {qqql of each generation) equals zero. The values of masses of the 4th 
generation particles we take from paper [ISj in which the fit to the electroweak observables 
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Figure 3: The final baryon asymmetry versus the initial asymmetry ub'/^ as a function 
of sphaleron freeze-out temperature Tgph for the unmixed fourth generation is shown by a 
dashed (blue) line. It is analogous to Fig. 2 from p3] but for m^r = 57.8 GeV, = 107.6 
GeV, rrit' = 634 GeV, mh> = 600 GeV. The final baryon asymmetry for the case of the mixed 
fourth generation quarks and the unmixed fourth generation leptons is shown by a solid 
(green) line. 



for higgs mass rriH = 125 GeV was performed and recent LHC bounds on the masses of t'- 
and 6'-quarks were taken into account: 

mt' = 634 GeV , nib' = 600 GeV , 
ruE = 107.6 GeV , mjv = 57.8 GeV . (38) 

The dashed blue line in Fig. 3 corresponds to the case of the unmixed fourth generation 
particles considered in [12] . The results for the case of the strongly mixed fourth generation 
quarks and the unmixed fourth generation leptons are shown by the solid green line. In order 
that leptons, A^, do not decay before the sphaleron freeze-out, which happens at tu ~ 10~^° s, 
the mixing angles of with three light neutrinos should be small: 6i < 10~^ (A^ decays 
through four fermion interaction). Assuming similar bound 6 < 10~^, the existence of heavy 
Dirac sequential neutrino with = 50 — 100 GeV is compatible with the search at LEP II 

According to the standard cosmological scenario nonrelativistic matter started to domi- 
nate the cosmic energy density at redshift z ~ 10^. If we demand that A^ should decay before 
that epoch, its life-time should be sufficiently short, tn < 10^^ s, from which we obtain the 
lower bound 9 > 10~^^. (Let us note that direct searches exclude A^ as a unique dark matter 
candidate [15]). 

A stronger bound on tn follows from the equilibrium form of the energy spectrum of 
CMB. According to ref. [E] a large influx of energy into the usual cosmological cosmic 
background would be thermalized if it took place before z ~ 10^. Otherwise the observed 
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black body spectrum of CMB would be noticeably distorted. Since the precision of the 
spectral shape is at the level of 10^^, only a very small distortion is permitted. 

The condition that decays before or at z ~ 10'' demands r^r < 10^ s, or > 10^^'^. 
If indeed decays before z ~ 10^, the contribution from its decay to the energy density 
of CMB would be not larger than 1% and the ratio of baryon to photon number densities 
r\B = UB/n^ at BEN epoch and at CMB recombination would be slightly different but in 
principle measurable by the light element abundances. 

More interesting and pronounced effect appears if heavy quarks of the 4th generation 
are long-lived. In this case we cannot explain the large value of CPV in D decays but may 
explain the difference of rjs at BBN epoch [tibbn) and at the recombination [rjrec) which 
is probably requested by the recent data on the light element abundances If heavy 

baryons of the 4th generation decays after BBN but before the hydrogen recombination, the 
number of light baryons in the comoving volume at BBN would be different from that at the 
recombination. The ratio riBBN/Vrec at these epochs could be either larger or smaller than 
unity depending upon the value of the baryon asymmetry in the heavy quark sector and the 
energy influx to CMB from the heavy baryon decays. So in principle both rise or decrease 
of riBBN is possibl^. 

In the limit Tsph — heavy particles of the fourth generation are not produced: = 
L4 = 0, -B' = L = 3A. In the physically interesting opposite limit Tgph ^ m^v the value 
of baryon asymmetry is nonzero since the right-handed neutrinos of three light generations 
are not produced in the primordial plasma violating symmetry between the leptons of four 
generations which would occur at T ^ rriN- The characteristic time of the right-handed 
neutrino to thermalize is T/mf, and for mj, < 1 keV (which is valid for three light neutrinos) 
this time is longer than the Universe age, tjj = Mpi/T^ for T = Tgph ~ 200 GeV, [11]. 

7 Conclusions 

In Introduction we determined what ratio of the penguin to the tree amplitudes of D ^ 
7r+7r~ decay is needed to get the observed CP asymmetry. In Section 2 we found that 
the factorization describes the tree amplitude with good accuracy; concerning the penguin 
amplitude it appears to be twenty times smaller than one needs to describe the experimental 
data on Acp- In Section 3 we demonstrated that in the case of -B — )■ n'^K^ decay the 
factorization underestimates the penguin amplitude by factor 2. In the case of Ks — )■ tt+tt" 
decay the penguin amplitude is enhanced by factor 2-3 in comparison with the factorization 
result. 

Thus if confirmed on larger statistics and future systematics result ([T]) demands New 
Physics. 

In Section 5 we demonstrated that the fourth quark-lepton generation may enhance the 
penguin amplitude describing the experimental data. If the leptons of the fourth genera- 
tion weakly mix with three light generation leptons, then the baryonic charge generated at 
high scale escapes the erasure by sphalerons and survives till now according to the results 
presented in Section 6. 

^ Since both the value of tjbbn and the number of Hght neutrino species influence nucleosynthesis, the 
change in the value of tibbn can be formulated as an additional (positive or negative) number of light 
neutrino species [IT]. 
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Below we derive equations used in Section 6 to find the dependence of the baryon asym- 
metry of the Universe on the sphaleron freeze-out temperature. In this Appendix we closely 
follow paper [T2] . 

Being interested in the values of the asymmetries at sphaleron freeze-out temperature 
we should assume that the electroweak phase transition already has occured and the neutral 
Higgs boson condenses. That is why the Higgs boson chemical potential is zero. Sometimes 
in the literature the baryon density in the electroweak unbroken phase is looked for. In this 
case the Higgs boson does not condense and its chemical potential is nonzero. To find it an 
additional equation is needed. It is provided by the condition that the density of charges with 
which the massless bosons interact should be zero, and in an unbroken phase there are two 
such charges: the hypercharge and the third projection of a weak isospin. The baryon density 
in the unbroken phase is analyzed, for example, in book [18] and it differs from its value in 
a broken phase. Since the right-handed components of quarks and leptons emitting neutral 
Higgs transform to the left-handed components the chemical potential of both components 
are equal: = /i^^ = fiu, fJ'd^ = l^d^ = f^d, /^e^ = /^e^ = /^e- The analogous relations are 
valid for the particles of the second and third families. The right-handed neutrinos of three 
light generations are not thermalized and should not be taken into account (see the end of 
Sect. 6). The fourth generation right-handed neutrinos, being heavy, rapidly thermalize: 
/iTVfl = fJ'NL = f^N- The chemical potentials of up and down weak isospin components are 
related by W~ chemical potential: fid = fJ'W + f^u, fJ^e = jJ'W + /^i^, I^e = l^w + I^^n- Mixing 
of quarks of four families and leptons of three families equilibrates the chemical potentials 
of the particles with the identical gauge quantum numbers. As a result four independent 
chemical potentials remain: /ztv, IJ^w and /i = n^,^ + /ij^j + jiy.^ = 3fi^. 

The particle number densities depend on their (Fermi or Bose) statistics, temperature, 
chemical potential, and masses. The chemical potential of an antiparticle is opposite to that 
of the particle. The asymmetries and, hence, chemical potentials are very small. Expanding 
the equilibrium integrals for the asymmetry over /j, we obtain: 



I (f ) '^/(^) ' if p is a fermion , 

where gp is the number of the degrees of freedom of the particle p {gq = gi = 2, g,^ = 1, 
gN = 2, gw = 3) and x = m/T. Functions a{x) are normalized in such a way that 
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ab{0) = Q;/(0) = 1. In what follows we take into account the nonzero masses of the particles 
of the fourth generation, of t-quark, and of ly-boson. 

The condition of electroneutrality of the primeval plasma looks as: 

2 1 
Q = 3-[2(a„ + + a* + 3 gPl^d + «s + ttf) + ttfe') 

X (/ivy + - 2[{ae + a^ + ar){nw + /^i.)] - 2aE{lJ,w + I^n) - 

- 3 ■ 2awf^w = , (A.2) 

(1 + 2at + 20:4/ — ay)iiu — (6 + ay + aE + 3aiy)/iiy — — cieI^n = 0. (A. 3) 

Here and below we omit irrelevant factor 

The sphaleron transition converts qqql combination of each generation into vacuum, which 
gives: 

12iiu + Sf^w + /i + /iTV = . (A.4) 

The remaining two equations are two superpositions of B', L, and L4 conserved under 
sphaleron transitions thus being equal to their initial values. The expressions for these 
quantities look like: 

L4 = 2aEfJ'E + 2aNjjN = 2{aE + aN)lJ'N + 2aEHw , (A-5) 
L = 2{ae + a^ + ar)fie + {c^u, + (Xu^ + «i/J^ = 3yu + 6fj,w , (A. 6) 

B' = 2 ■ 3 ■ ^[{au + ac + at + at')fiu + {<yd + as + C(b + C(b')lJ'd\ = 

= 2{2 + at + at') fiu + 2{3 + ab') {fin + fiw) ■ (A.7) 



Thus we have four equations which determine the chemical potentials: ( ]A.3I) . ( ]A.4I) . and 
the remaining two: 

B' — L — L/^ = 2(5 + at + a^ + ab')fJ.u + 2{ab' — aE)liw — 

—3ji — 2{aE + aN)nN = ^ , (A-8) 
L — 3L4 = 6(1 — a;£;)/ivy + 3/i — 6(a£; + a;Ar)/iAr = 3A , (A. 9) 

where we take the initial values analogous to those of ref. [12]: B'^ = Lq = 3 A and L4 = 0. 
When temperature is much larger than the masses of all the particles, all Q!j are equal to 
one we obtain: 

11 



B' 

'a 



(A.IO) 
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If the right-handed neutrinos of three light generations thermalized then the equation 
flA.6l) would be substituted by 

L = Afi + 6fiw , (A.ll) 
and the baryon asymmetry at T ^ would vanish. 
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